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ON THE STABILITY OF SELF-SIMILAR SOLUTIONS TO 
NONLINEAR WAVE EQUATIONS 

OVIDIU COSTIN, ROLAND DONNINGER, IRFAN GLOGIC, AND MIN HUANG 


Abstract. We consider an explicit self-similar solution to an energy- 
supercritical Yang-Mills equation and prove its mode stability. Based 
on earlier work by one of the authors, we obtain a fully rigorous proof of 
the nonlinear stability of the self-similar blowup profile. This is a large- 
data result for a supercritical wave equation. Our method is broadly 
applicable and provides a general approach to stability problems related 
to self-similar solutions of nonlinear wave equations. 


1. Introduction 

The development of singularities in finite time is one of the most stun¬ 
ning features of nonlinear evolution equations. Singularity formation (or 
“blowup”) of the solution signifies a dramatic change in the behavior of the 
underlying model or even the complete breakdown of the mathematical de¬ 
scription. On the level of a fundamental physical theory, blowup occurs in 
Einstein’s equation of general relativity to indicate the dynamical formation 
of a black hole. However, a rigorous treatment of Einstein’s equation in this 
context is hopeless at the present stage of research. Consequently, it is a 
reasonable strategy to resort to simpler toy models that capture some of the 
features of the more complicated system. Natural candidates in this respect 
are energy-supercritical nonlinear wave equations with a geometric origin 
such as wave maps or Yang-Mills models. 

The easiest way to demonstrate finite-time blowup in a given evolution 
equation is to construct self-similar solutions. In exceptional cases it is even 
possible to obtain closed-form expressions. The relevance of such solutions 
depends on their stability. After all, one would like to obtain information 
on the generic behavior of the system. However, already at the linear level 
the stability analysis of self-similar solutions to nonlinear wave equations is 
very challenging since one is confronted with highly nonself-adjoint spectral 
problems. Consequently, standard methods do not apply. This fact poses a 
serious obstacle to any rigorous analysis of the blowup dynamics. 
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In the present paper we develop a general approach which is capable 
of handling the difficult nonself-adjoint spectral problems related to self¬ 
similar blowup. For the sake of simplicity, however, we focus on the concrete 
example of an energy-supercritical Yang-Mills equation that displays blowup 
via an explicitly known self-similar solution. 


1.1. An energy-supercritical Yang-Mills model. For /i G {0,1, 2,... , 5} 

let Afj_ : —?■ so(5) be a collection of five fields on (1 -|- 5)-dimensional 

Minkowski space with values in the matrix Lie algebra of SO(5). In other 
words, for fixed /x and {t, x) G x) is a skew-symmetric real (5 x 5)- 

matrix. One sets 

Ffiu dfj^Ap duA^ -\- [A^, Ay\ 

and considers the action functional^ 

[ tv{F,,Fn- ( 1 . 1 ) 

jRl.S 

Formally, this is reminiscent of Maxwell’s theory. However, the commuta¬ 
tor in the dehnition of F^i^ introduces a very natural nonlinearity. In this 
sense, Yang-Mills theory can be viewed as a nonlinear generalization of elec¬ 
trodynamics. The Euler-Lagrange equations associated to the action (1.1) 
are 

+ =0 

and the ansatz [17, 5] 

AF{t,x) = {6f,x^-6lx^)^/^^^ 

yields the scalar nonlinear wave equation 


2^,. , 3V’(i/’ + l)(iA + 2) 

Wtt Yrr “1“ 9 — 0? 


( 1 . 2 ) 


'ip = for the auxiliary function i/; : M x [0, oo) —>■ M. Eq. (1.2) has 

been proposed as a model for singularity formation in Einstein’s equation 
[4, 5, 2, 21]. In general, (classical) Yang-Mills fields attracted a lot of interest 
by both the physics and mathematics communities, see e.g. [1, 18, 19, 22, 
23, 3, 30, 25, 29, 11, 24, 27, 28]. 

Eq. (1.2) is energy-supercritical [4] and large-data solutions can develop 
singularities in finite time as is evidenced by the existence of self-similar 
solutions of the form = /(j^), see [7]. Bizoii [5] found an explicit 

example of this kind given by 


Mt, r) = /o(t^), Mp) = 

Numerical investigations [4, 5, 2] yield strong evidence that the solution ipQ 
gives rise to a stable self-similar blowup mechanism. Motivated by this. 


^Einstein’s summation convention is in force. Greek indices take the values 0 to 5 
whereas latin indices run from 1 to 5. Our convention for the Minkowski metric is 


7? = diag(-l, 1,1,1,1,1). 
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the second author [13] developed a complete nonlinear stability theory for 
the solution ipQ, see also [12, 16, 14, 15] for other types of nonlinear wave 
equations. However, the results in [13] are conditional in the sense that they 
depend on a spectral assumption which could not be verified rigorously so 
far. It is the aim of the present paper to close this gap. 


1.2. The mode stability problem. The first important step in a stability 
analysis of the solution -00 is to rule out unstable modes. To this end, one 
introduces similarity coordinates [2] 

r = -log(l-t), p=j^. 

Eq. (1.2) transforms into 

(pTT “1“ “t“ ^pYrp (1 P )\Srpp “1“ p^p) ^2 — ^ (^‘^) 

where 4'{t,p) = p(l — e~'^,e~'^p). Due to finite speed of propagation one 
is mainly interested in the behavior inside the backward lightcone of the 
singularity, which corresponds to the coordinate domain r > 0, p € [0,1]. 
Note that the self-similar solution is independent of r and simply given by 
fo{p)- Next, one inserts the mode ansatz 

0(d p) = fo{p) + e^'^uxip), A € C 


and linearizes in u\. This yields the ODE spectral problem 

— (1 — p^){u'^ + ^u'x) + 2Xpu'x + A(A -|- l)ttA H —~ ^ 
for the function ua, where the potential V is given by 


(1.4) 


H(p) = 6 + 18 /o(p) + 9/o(/9)2 = 6 


25 - 90p^ + 33/9^ 

(5 + 3p2)2 


Observe that Eq. (1.4) has a singular point at the lightcone p = 1 which is 
a consequence of the fact that lightcones are the characteristic surfaces of 
Eq. (1.2). 

Admissible solutions of Eq. (1.4) with Re A > 0 lead to instabilities of /o 
at the linear level. However, it is not entirely trivial to determine what “ad¬ 
missible” in this context means. This question can in fact only be answered 
once one has a suitable well-posedness theory for Eq. (1.3). The necessary 
framework is developed in [13] and it turns out that if Re A > 0, only smooth 
solutions are admissible. Consequently, a nonzero solution u\ € C°°[0,1] of 
Eq. (1.4) with Re A > 0 is called an unstable mode. The corresponding A is 
called an (unstable) eigenvalue. As a matter of fact, there exists an unstable 
mode. The function 

Mp) := 


turns out to be a smooth solution of Eq. (1.4) with A = 1, as one easily 
checks. However, this mode is not a “real” instability of the solution /o but 
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rather a consequence of the time translation symmetry of Eq. (1.2). Indeed, 
the prohle /o defines in fact a one-parameter family of blowup solutions 
given by 

where T > 0 is a free parameter. By the chain rnle it follows that 

solves the linearized equation. These observations lead to the following 
definition. 

Definition 1.1. The solution tpQ (or fo) is said to be mode stable if ui is 
the only unstable mode. 

1.3. The main result. With these preparations at hand we can formulate 
our main result. 

Theorem 1.2. The self-similar solution tpQ is mode stable. 

The hrst result of this kind was proved very recently for a similar problem 
related to the wave maps equation [10]. However, the method we develop 
here is different and much more effective. As a consequence, the main ar¬ 
gument fits on a few pages and the method easily generalizes to other types 
of nonlinear wave equations. In view of the fact that rigorous research on 
self-similar blowup in supercritical wave equations was blocked for a long 
time by the difficulties related to these spectral problems, we hope that our 
method will trigger new developments in the field. In this respect we also 
remark that Theorem 1.2 in conjunction with the theory developed in [13] 
yields a fully rigorous proof of stable self-similar blowup dynamics for the 
Yang-Mills equation (1.2). The precise statement is given in [13], Theorem 
1.3. We emphasize that this is a large-data result for an energy-supercritical 
wave equation. 


2. Removal of the symmetry mode 

Although the eigenvalue A = 1 is not connected to a real instability of the 
solution fjQ, it is still inconvenient for the further analysis. Consequently, it 
is desirable to “remove” it. This can be done by a suitable adaptation of a 
well-known procedure from supersymmetric quantum mechanics which we 
recall here briefly. 

2.1. Interlude on SUSY quantum mechanics. Consider the Schrddinger 
operator H = —8“^ -|- U on L^(M) with some nice potential V and suppose 
there exists a ground state fo € L^(R) PI C'°°(M), i.e., /q = Vfo. Assume 
further that fo has no zeros. Then one has the factorization 
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By interchanging the order of this factorization, one defines the SUSY part¬ 
ner H of H, i.e., H := QQ*. Explicitly, the SUSY partner is given by 

=---dl+y 

fl2 

where V = —V -|- 2^ is called the SUSY potential. The point of all this 

JQ 

is the following. Suppose A is an eigenvalue of H, i.e., Hf = Q*Qf = Xf 
for some (nontrivial) /. Applying Q to this equation yields QQ*Qf = XQf, 
i.e., HQf = XQf. Thus, if Qf / 0, i.e., if / ^ kerQ, A is an eigenvalue of H 
as well. Obviously, we have ker Q = (/o) and thus, if A / 0 is an eigenvalue 
of H, then it is also an eigenvalue of H. Moreover, 0 is not an eigenvalue of 
H for if this were the case, we would have QQ* f = 0 for a nontrivial /, i.e., 
/ € kerQ* or Q*f G kerQ. The former is impossible since kerQ* = (^) 
but ^ L^(M). The latter is impossible since TgQ* T kerQ. In summary, 

H has the same set of eigenvalues as H except for 0. 


2.2. The supersymmetric problem. Now we implement a version of this 
SUSY factorization trick for our problem. Note that the Frobenius indices 
of Eq. (1.4) at yO = 0 are {—3,2} and at /? = 1 we have {0,1 — A}. Suppose 
ux is an unstable mode of Eq. (1.4) and A 7 ^ 0. By definition, ux G C'°°[0,1] 
and from Frobenius theory it follows that \ux{p)\ ~ as p —)• 0-|- as well as 
|u;v(/ 9 )| ~ 1 as p —>■ 1—. We define a new function vx by^ 

ux{p) = P~\^ - P^)~^^‘^vx{p). 


From Eq. (1.4) it follows that vx satisfies 

// y{p) 

- < + 2 n = 

p^{l- pQ 

For A = 1 we have 

vi{p) = p{l - p‘^)^ui{p) = 
We rewrite Eq. (2.1) as 


A (2 - A) 

80p3(l -p2)l 
(5 + 3p2)2 


-v'l + Vivx 


A (2 - A) - 1 


Vx 


( 2 . 1 ) 


^Observe that this transformation depends on A. This is the reason why Eq. (1.4) is 
not equivalent to a standard self-adjoint Sturm-Liouville problem. What happens is the 
following. Since |ma(p)| — 1 as p —>■ 1 — , the corresponding v\ behaves like |nA(p)| — 
(1 — The Hilbert space in which the spectral problem for va is symmetric is 

LJ,(0, 1) with the weight iu(p) — ^ 7375 ^. Thus, if Re A < 1, the admissible solution vx 
does not belong to L^(0, 1)! Consequently, for Re A < 1 the self-adjoint formulation does 
not yield any information. This shows that the spectral problem (1.4) is truly nonself- 
adjoint in nature. In particular, there can be nonreal eigenvalues. For Re A > 1, on the 
other hand, one can indeed use Sturm oscillation theory to exclude eigenvalues. 
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with 

.. r ^ yjp) 1 

/02(l-p2) (l-p2)2' 

Then we have v'l = Vivi and thus, Eq. (2.1) may be factorized as 


{-dp- 




A (2 - A) - 1 
(1 - 


v\ 


or 

-(1 - p^f{d, + = [A (2 - A) - 

We set v\ = {dp — ^)fA and apply the operator ^ to the equation 
which yields the supersymmetric problem 


- {dp - J)[(l - p^f{dp + $)]hA = [A (2 - A) - IJuA. ( 2 . 2 ) 
Note the asymptotics 

±{p) = 2,p-^ + 0{p) (p^0+) 

Vl 

^(p) ~ -d(l-p)“^ (P^l-). 

Consequently, from the representation v\{p) = p^h\{p^), where h\ is ana¬ 
lytic near 0, we get v\{p) = 0{p^) near p = 0 and from v\{p) ~ c (l-p)V2 
we infer v\{p) ~ c {1 - near p = 1 (unless A = 1). Writing out 

Eq. (2.2) explicitly yields 

- (1 - Pfvl + 4p(l - P)iii + = A(2 - \)v, (2.3) 

with the supersymmetric potential 


V{p) = 20 


15 - 2p2 3p^ 

(5 + 3p2)2 


Setting u\{p) = p ^(1 — p^Y ^^^w(p) we find the equation 


— (1 — p'^){u'^ + ^u'x) + 2\pu'x + (A^ -|- A — 2)12^ H-— 0. (2.4) 

P p^ 

Note that the Frobenius indices of Eq. (2.4) are {—4, 3} at 0 and {0,1 — A} 
at p = 1. With minor modifications the same procedure can be performed 
in the case A = 0. As before, we say that A € C is an unstable eigenvalue of 
Eq. (2.4) if Re A > 0 and there exists a nontrivial solution u\ € C'°°[0,1] of 
Eq. (2.4). In summary, we have proved the following result. 


Proposition 2.1. Let X Y ^ be an unstable eigenvalue of Eq. (1.4). Then 
A is an unstable eigenvalue of Eq. (2.4). 
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3. Absence of unstable eigenvalues for the supersymmetric 

PROBLEM 

In this section we exclude unstable eigenvalues of Eq. (2.4). Via Proposition 
2.1 this implies the main result Theorem 1.2. 

Theorem 3.1. The supersymmetric problem Eq. (2.4) does not have unsta¬ 
ble eigenvalues. 

The Frobenius indices of (2.4) at 0 are —4 and 3, hence the solution 
analytic at 0 has the power series representation 

CX) 

ao^O. (3.1) 

n=0 

Note that A is an eigenvalue of (2.4) if and only if the radius of convergence 
of (3.1) is greater than 1. Therefore, our aim is to prove that for any A 

in the closed right half-plane (which from now on we denote by H), (3.1) 

cannot be analytically extended through p = 1. 

By substituting (3.1) into (2.4) we obtain a four term recurrence relation 
(with the initial condition uq = 1 and = 0 for n < 0) 

P3{n)an+3 + P2{n)an+2 + Pi{n)an+i +po{n)an = 0, (3.2) 

where 

P 3 (n) = -lOOn^ - 950n - 1950, 

P 2 (n) = -2071^ (lOOA - 150)n -h 25X^ + 375A - 370, 

pi{n) = 84/1^ -h (120A + 462)n -h 30A^ -h 330A + 630, 

Po{n) = 36n2 + (36A + 126)n -h OA^ + 63A + 90. 

One can check that a„ = (—3/5)” is an exact solution to (3.2), hence the 
order of the recurrence (3.2) can be reduced by one through the substitution 

bn — O.n+1 T g an. (3'3) 

This yields a three term recurrence relation for bn 

q 2 {n)bn +2 + qi{n)bn+i + qo{n)bn = 0, (3.4) 

where 

q2{n) =P3{n), 

qi{n) =P 2 {n) - |p3(n), 

qo{n) =Pi{n) - f P2(n) + ^Psin). 

After substituting for pi{n) in the last three relations, dividing all of them 
by 5 and using the qi notation for the new coefficients, we get 

q 2 {n) = —20n^ — 190n — 390, 

qi{n) = -k (20A -k 84)n -k hA^ 75A -k 160, 
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qo{n) = 12n^ + (12A + 42)n + 3A^ + 21A + 30. 

By letting = qi{n)/q 2 {n) and Bn = qo{n)/q 2 {n), (3-4) becomes equivalent 
to 

bn+2 4 “ -^nbn+l 4 “ Bnbn — 0 , ( 3 . 5 ) 

with the initial condition b -2 = 0 and 6_i = 1. 


Lemma 3.2. Given A in the complex plane, either 


lim 

n—^oo 


b-n+l ('^) 

“MAT 


= 1 , 


or 


lim 

n—j-oo 


^n+1 ( A) 

~b;w 


3 

5’ 


(3.6) 

(3.7) 


Proof. Since lim„^ooA„(A) = —2/5 and lim„^ooB„(A) = —3/5, the char¬ 
acteristic equation associated to (3.5) is 

It-1 = 0. (3.8) 

As the solutions to (3.8) (1 and —3/5) have distinct moduli, by a theorem of 
Poincare (see, for example, [20], p. 343, or [6]), either bn is zero eventually 
in n, or lim„_).oo &n+i(A)/&n(A) exists and it is equal to either 1 or —3/5. 
Now, for a fixed A, bn cannot be zero eventually in n, since by backward 
induction from (3.5) one would get b-i = 0, hence the claim follows. □ 


Note that in order to prove Theorem 3.1, it suffices to show that (3.6) 
holds for all A in H, for that implies non-analyticity of (3.1) at 1. Indeed, 
defining fx by (3.1) and gx by gx{p) = ao(A)p + X]^o ^ri(A)p^’^+^, one easily 
checks that 

f\{p) = ^p2\^ 9x{p)- (3.9) 

So if (3.6) holds and therefore gx is singular at 1, then, by (3.9), so is fx. 
Let rn = bn+i/bn. Then from (3.5) we obtain 


I’n+l — An 



(3.10) 


where 


r_i = ^ = -A_ 2 (A) = —A2 + ^A + —. 
6_i ^ ^ 18 18 15 


(3.11) 


The idea is to find a “simple”, provably close approximation to rn in H, that 
converges to 1 for any fixed A, which would then imply (3.6). 

We use the quasi-solution approach, initially developed for ordinary dif¬ 
ferential equations in [8, 9], which we here, in a sense, extend to difference 
equations of type (3.10). Namely, as a quasi-solution to (3.10) we define 


~ lAl ^ n 4“ 2 

”■ 4n^ 4- 31n -|-43~*~n4-4~*~n4-4 


(3.12) 
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Of course, the choice is not arbitrary, and in §4.1 we describe in some detail 
how to obtain such an approximate solution. The quasi-solution turns 
out to be a good approximation to in the whole of El. 

Lemma 3.3. ri and for n > 1, are analytic in H. 

Proof. From (3.10) and (3.11) we compute 

1 25A® -t 825A5 -h lOOlSA^ -F 69735A3 -F 207694A2 -F 260856A -F 96192 
’~ 78 25A4 -F 450A3 -F 2735A2 -F 5070A -F 2016 ’ 

The denominator of ri and the polynomials rn(A) for n > 1 are Hurwitz- 
stable i.e., all of their zeros are in the (open) left half-plane, which can 
be straightforwardly checked by, say, the Routh-Hurwitz criterion or its 
reformulation by Wall (see [32] or §4.2)^. The conclusion follows. □ 


Now, let 

(3.13) 

Substitution of (3.13) into (3.10) leads to the following recurrence relation 
for 6n, 

dn+l = + Cn- — ' —( 3 - 14 :) 

1 + C>n 

where 


.^nrn Bji A 

e„ =-1 and = - 


rnrn+i 




(3.15) 


Lemma 3.4. The following estimates hold in El, 



1 

l^nl < 

l<^n| < i 

5 


n > 1, 
n > 1. 


(3.16) 

(3.17) 

(3.18) 


Proof. The method of proof is the same for all three quantities, so we illus¬ 
trate it only on Cn. 

Lemma 3.3 and (3.15) imply that Cn is analytic in El. Also, being a ratio¬ 
nal function, Cn is evidently polynomially bounded in El. Hence, according 
to the PhragmOT-Lindelof principle”^, it suffices to prove that (3.18) holds on 
the imaginary line. To that end, we first bring C'„+i(A) to the form of the 
ratio of two polynomials Pi{n, A) and P 2 {n, A)^. Then, for t real, |C,i+i(it)p 
is equal to the quotient of two polynomials, Qi{n,t‘^) = |Pi(n,it)p and 
Q 2 in,t‘^) = \P 2 {n,it)\‘^. In order to show that |C'„+i(it)| < 3/5, for all real 


^There are, of course, elementary ways of proving this claim. However, the suggested 
approach is more general. 

^We use the sectorial formulation of this principle, see, for example, [31], p. 177. 

^For all three quantities, straightforward calculations would lead to the form that we used. 
However, to prevent possible ambiguity, in §4.3 we give the explicit form (as a ratio of 
polynomials) for all three quantities. 
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t and n > 0, all we need is to show that \Cn+i{it)\‘^ = Qi{n,t‘^)/Q 2 {n,t'^) < 
9/25, or equivalently 9/25 ■ Q 2 — Qi >0. Using elementary calculations, we 
see that 9/25 ■ Q 2 — Qi has manifestly positive coefficients, and the variable 
t appears with even powers only. Thus, (3.18) holds on the whole imaginary 
line, and the result follows. □ 

Proof of the Theorem 3.1. Prom (3.14) and Lemma 3.4, a simple inductive 
argument implies that 


|<5n| < for all n>l, and A € HI. (3.19) 

Since for any fixed A, lim„_).oo rn(A) = 1, (3.13) and (3.19) exclude the 
possibility of (3.7). Hence, (3.6) holds in H, and the claim follows. □ 


4. Appendix 

4.1. Description of how to obtain a quasi-solution. First, the minimax 
polynomial approximation® of degree two to over an interval [0,10] is 
found, where n ranges from 0 to 20. Then, appropriate rational functions 
in n are fitted to the coefficients of the approximation polynomials. 

We should point out that interval of polynomial approximation and the 
range of values of n can vary, and the ones from the description are just 
our choice. We choose quadratic polynomial approximations due to the fact 
that is a ratio of two polynomials whose degrees differ by two. 

4.2. Wall’s criterion for Hurwitz-stability. Let P{z) = z^ + aiz'^-^ + 

■■■ + On he a polynomial with real coefficients, and let Q{z) = + 

a^z'^~'^ + • • • be the polynomial that contains exactly those terms of P{z) 
that have odd-indexed coefficient. Then all the zeros of P{z) have negative 
real parts if and only if the quotient Q{z)/P{z) can be represented in a finite 
continued fraction form 


1/(oi -|- 1/(02 + 1/(03 + ... + 1/ On) . . .), 

where oi = ciz + 1, 02 = C 2 Z, ..., an = CnZ, and the coefficients ci, C2, ..., 
Cn are all positive. 

In our case, for the denominator of ri, the coefficients Cj are ci = 1/18, 
C2 = 135/736, C3 = 33856/64863 and C4 = 36035/15456, and for f„, ci = 
(n + 4)/(4n^ -t 31n -t 43), and C2 = l/(n -t 2). 


®The minimax polynomial approximation of degree n to a continuous function / on a 
given finite interval [a, b] is defined to be the best approximation, among the polynomials 
of degree n, to / in the uniform sense on [a, b]. For the proof of existence and uniqueness 
of this approximation and an algorithm to obtain it, see [26], §2.4. 
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4.3. Detailed expressions for Cn, £n and (ii. We give details of these 
quantities in order to fully clarify the notations. We have 

Cn+i = Pi{n, X)/P 2 {n, A), 

where 

Pi(n, A) = - 3(n + 5)(n + 6){W + 39n + 78)(dn^ + 47n + 121) 

X [A^ + (4n + 11)A + dn^ + 22n + 28] 

and 

P 2 {n, A) =10(2/1^ 23n + 60)[(n + 5)A2 + (dn^ + 39n + 78)(A + n + 3)] 

X [(n + 6 )A^ + (4n^ -|- 47n + 121)(A + n + 4)], 
respectively. Furthermore, e„+i = P 3 {n, X)/P 2 {n, X), where 

Ps^n, A) =5(n + l)(n + 5)(n -|- 6)A'^ 

- 5(8n'^ + 158n^ + 1095n2 -P 3171n + 3162)A3 

- (112n^ + 2364n^ 17243n^ + 4880571^ + 33244n - 36060)A^ 

- 4(4n^ 39n + 78) (4n^ + 47n + 121) 

X [(3n^ + 5n - 3)A - dn^ - 3n + 36]. 

Finally, 

-5A2(15A3 - 20A2 - 939A + 1412) - 36(1093A - 256) 

^ “ (5A2 78A + 234)(25A4 450A3 2735A2 5070A + 2016)' 
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